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ABSTRACT 

. - . The Lamb shift measurement and theory are now both a dynamicaUy developing field 

' and we give a review of the current data. Critical comparison of theory and experiment can 

, be done using a value of the proton charge radius and we pay attention to several results 

' of its determination. 

< 

The talk is devoted to the Lamb shift in the hydrogen atom. The workshop is for 
exotic atoms and first of all I would like to say that to my mind the hydrogen is one of 
them. The possibility both to calculate and to measure different energy intervals with 
^ I extremely high accuracy make hydrogen a quite exotic system. Next , one has to remember 
' that there is actually no special separate theory of the positronium or muonium atom, 
^ ■ investigations of which were discussed here in detail. The same expression can be of use 
^ . in a number of calculations for several atomic systems (see e. g. Refs. i. One more 
Q \ point which is quite similar to the pionium and pionic hydrogen is that we investigate an 
O I atomic system but after all the result is important for particle physics, namely for the 
^ ' proton charge radius. All of these reasons show that hydrogen atomic properties are to 

be discussed among exotic ones. The workshop talk is based on a Max-Planck-Institut 
Oh. fiir Quantenoptik report i and all references can be found there. 

r-| ! The hydrogen atom is one of the most important QED systems. In contrast to 

muonium and positronium the nucleus is a proton and hence the energy of the level is 
influenced by the proton structure, i. e. by the strong interaction. The knowledge of 
the proton radius leads now to a limit of the theoretical value for the Lamb shift. We 
discuss here the problem of the radius determination, the most popular values of which 
are summarized in Table 1. We give a critical review of all of these values, as well as of 
the theory and the experiment on the Lamb shift. 

Let us briefly describe all items of the Table 1. The elastic electron-proton scat- 
tering data were obtained with small momentum transfer and were extrapolated to zero 
momentum. The other way to extract the radius from such data is based on a dispersion 
relation approach which allows to involve data from other kinematic areas into fitting. 
The numerical calculation within the chiral limit of the lattice QCD can be corrected due 
to the chiral perturhative theory. And indeed one of the ways to determine the radius is 
based on the hydrogen Lamb shift investigation. First we mention that we expect that the 
uncertainties presented in Table 1 accordingly to the original works are significant 

underestimations and we will not consider those anymore here. The details can be found 



X 
J3 



m our review 



'•'E-mail: sek@mpq.mpg.de; sgk@onti.vniim.spb.su 



1 



Value 


Reference 


Method 


0.809(11) fm 
0.862(12) fm 
0.64(8) fm 
0.88(3) fm 
0.847(9) fm 
0.890(14) fm 


Stanford, 1963 i 
Mainz, 1980 1 
Draper et ai., 1990 1 
Leinweber, Cohen, 1993 
Mainz, 1996 i 
Garching, 1997 1 


scattering & empirical fitting 
scattering & empirical fitting 
lattice QCD in chiral limit 
lattice QCD & chiral perturbation 

dispersion relation fitting 
hydrogen Lamb shift measurements 



Table 1: Proton charge radius 



We start our consideration with the Lamb shift. Several results obtained within 
different approaches are summarized in Table 2 (see Fig. 1 for more detail). The results 
presented are found within different experiments: 

• the LS result is an average value of the best direct measurement of the Lamb splitting 
S by Lundeen and Pipkin 113 and older works. We also take into account a recent 
Daper by van Wijngaaden, Holuj and Drake 0, published after our work on review 
13 was completed; 

• the LS/T result of an indirect measurement of the splitting S. The measurement 
was performed by Sokolov and Yakovlev 111 for the ratio of S and the 2pi/2 radiative 
width. The width was recalculated by Karshenboim recently including a leading 
radiative correction of relative order a{Zay\n{Za) (see Ref. Ill for detail) 

r(2pi/2) = '^a'Ry^i^l + \nQ {Zaf + a{Zaf (1.5084...)^ (^} ' 

• in case of FS the measured value is the fine structure splitting 2si/2 — 2p3/2, the 
most recent result for which was found by Hagley and Pipkin IIJ; 

• the OBF abbreviation means optical beat frequency. The value is an average one 
from recent Garching 111, Yale HI and Paris O data; 

• the CAF result is obtained by comparison of two optical frequencies measured 
separately, namely the Is — 2s interval from Garching H and 2s — 8s/ d from Paris 

SB. 

To recalculate the data obtained by FS, OBF and GAF methods one has to use a 
significant piece of theory. In case of the fine structure that is the theory of the 2p states 
and in case of the optical measurement that is theory of a specific difference 

A{n) = El{Is) - Eiins) . (1) 

First we explain the importance of this difference and next we consider the status of 
its calculation as well as of Ei{2pj). The progress in measurement of the either the 
Lamb splitting or the fine structure has been relatively slow. In contrast to that great 
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Method 


Lamb splitting 


LS 


1057850(7) kHz 


Lb/V 


1057858(2) kHz 


FS 


1057840(11) kHz 


OBF 


1057843(7) kHz 


CAF 


1057853(4) kHz 



Table 2: Experimental result for the Lamb splitting S = E{2si/2 — 2pi/2) 

development was obtained in optical measurement of the transition frequency between 
levels with different value of the principal quantum number n {gross structure). The 
highest precision was reached in two-photon Doppler-free transitions like Is — > 2s. The 
problem of utilizing those results was due to the Rydberg constant determination. There is 
no way to find it except by investigation of the gross structure. Hence to find anything for 
the Lamb shift one has to measure two optical transitions and to construct some difference 
in which the Rydberg contribution is canceled. It is possible to instrumentally extract a 
beat frequency directly or indirectly within an experiment. It is also possible to do that 
with data obtained in two independent determinations of different intervals. But still one 
has to solve one more problem: the Lamb shift result is after all a combination of the Lamb 
shift of the Is and 2s states and also a portion of a higher excited levels contribution is 
included. In order to manage that a specific difference A{n) was introduced by us some 
time ago (see Ref. for detail). This difference as well as the 2p state Lamb shift has 
a much better theoretical status than the ground state Lamb shift (or S). A number of 
contributions like e. g. a three-loop term which have not been known up to date for S 
are known for the difference and for the 2p states. 

The theoretical expression for the difference of Eq. (|l|) is of the form 

^ ' TT m2 \ 3 ko{ns) \ Mj 

(^4( ln(n) - ^(n + 1) + z^(2)) - ^^^^;^) In ^ + A^J'in) + G'J{Za) 
lAZm ( , , n-l\l a^iZafm , o 1 

-TIT + ^) - - '"("' + ^) h Tzaf ■ 

The results of the 2p states are similar to that 



+ {Zaf X 



3 



0.847 
0.862 

0.805 



Theo 



P1993 

Paris 
Garching 
Yale 

HP 



P1997 



CAF 



OBF 



FS 



Old 



WHD 
LP 



LS 



SY / K ( Ls/r ; 

Old ^ 



1057800 Lamb splitting S [ kHz ] 1057900 



Figure 1: Different values of the Lamb splitting S 
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The one- loop contribution due to the vacuum polarization {Aq^^^) was found for an arbi- 
trary ns state by Karshenboim and Ivanov and for an npj state by Manakov, Nekipelov 
and Feinstein. The self-energy term {G^^{Za)) is mainly determined by a value of Aqq^\ 
The last was found by Pachucki for Is and 2s and by Jentschura and Pachucki for 2p. 
The result has to be found extrapolating numerical data calculated by Mohr and by Kim 
and Mohr for higher nuclear charge Z. It is simpler to extrapolate the data for the differ- 
ence of Eq. and for the p state in comparison to extrapolation for the ground state. 
The leading two-loop logarithmic term was found by Karshenboim. Pachucki and Grotch 
found that the recoil term of the order [Zaf'w? jM for an ns state is scaling by Xjr? 
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and so it cannot contribute to the difference. In case of the 2p states this correction was 
evaluated by Golosov, Yelkhovsky, Milstein and Khriplovich. 

The theoretical progress mentioned for the A(n) and the 2p Lamb shift allows to 
obtain results in Table 2 with pretty small theoretical uncertainty . The QED results for 
4p states are also needed for the evaluation of some experimental data. The value of 
Ago was found by Jentschura, Mohr and Soff and the recoil effects were investigated 
by Yelkhovsky and Pachucki independently. That is, however, only a not significant part 
of the QED calculations which have to be done on a way to reach a value of the proton 
radius from hydrogen atom spectroscopy. It is necessary to calculate the Lamb shift of 
the 2s state (as we can see from the discussion above the significant part of the QED 
theory of Is and 2s states is the same and we speak here about the ground state). 

The expression for the ground state Lamb shift is much more complicated. The 
term most important for the further discussion is of the well known form 

dEnud{ns) = . (3) 

That includes the proton charge radius and subtracting the QED result (Eq. excl.) 
from the experimental value we can determine the proton charge radius. A detailed 
review of the ground state Lamb shift was presented by Pachucki et al. 0. Here we 
mention some important features. The most recent corrections considered there are the 
one-loop self-energy {Pachucki; Mohr), the two-loop contribution of the order a^(Za)^m 
{Pachucki and Eides, Grotch and Shelyuto), the leading two-loop logarithm of a higher 
order {Karshenboim) , a pure recoil term {Fell et ai; Grotch and Pachucki and Shabaev 
and coworkers) and the radiative-recoil correction {Bhatt and Grotch contradicting to 
Pachucki). 

We mainly agree with a consideration in review 111. However, our result is shifted 
by -3.6 kHz because the (ZQ;)^m log^ Za-correction has not been included there. We 
would also like to present here all sources of the theoretical uncertainty for the Lamb 
splitting S: 

• unknown a{Zaym and higher order one-loop corrections are estimated to 1 kHz; 

• a^{Za)^m\o^ Za and higher order two-loop terms can give up to 2 kHz; 

• the three-loop a^{Za)^m contribution is here estimated preliminary to 2 kHz and 
needs more understanding. 

The one-loop term is going to be calculated exactly and so the uncertainty is being 
lool 

removed 1^. 

Now let us turn to the discussion of the scattering results. We start from simple 
estimates. The value extracted from elastic electron-proton scattering cross sections is 
the electric form factor of the proton G{q^). We need to investigate it at low momen- 
tum transfer where the 'signal' which is G — 1 is mainly determined bv the radius term 
{{q'^Rp)/6). The G — 1 term lies at the condition of experiment in Ref. H between 1% and 
15% (see Fig. 2). The scattering radius in Table 1 is claimed i to be with uncertainty 
within about 1%. That means that all shifts on the level of a few percents are impor- 
tant. We discussed in our review i some QED corrections to the cross section which are 
expected to be on a few percent level and which are beyond evaluation of Ref. I. In 
this paper we concentrate our attention to more important problem. That is due to the 
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Figure 2: Experimental data for the proton electric form factor 

normalization of the data. It is clear that the value of the electric form factor at zero 
momentum transfer is equal to one 

Gth{0) = l. 

That is absolutely correct theoretically , but the form factor cannot be measured straight- 
forwardly. It is possible to measure only some cross section. Evaluating the data one can 
extract form factor from the Rosenbluth formulae. On this way the extracted experimen- 
tal value which is expected to be the form factor is only consistent with the true form 
factor within some uncertainty . As one of results the value of Gexp{0) is consistent with 
one and it has to be found within some fitting procedure. In other words we have to write 
a trivial equation 

where the value of Oq is consistent with one, but not equal to that a priori. For application 
to the low momentum transfer electron-proton scattering H we can expect that Oq is a 
constant (but it realy depends on experimental setup). 

In both Mainz papers i' i some special a priori prescriptions for were used. 
The problem of normalization was investigated by Wong il prior to publication of the 
dispersion paper i. He found that the result strongly depends on an assumption on 
the value of the constant Oq. The correct value of the radius corresponds to the free 
normalization (i. e. Oq is one of the fitting parameters). The constant 

ao = 1.0028(22) 
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Figure 3: Some fitting of tlie Mainz data witli different normalization constant 

found from tfie Mainz scattering data i a 'posteriori IH (see Fig. 3) is consistent witli one 
(tlie grescription of Ref. i is ao = 1) and witli the value of ao = 1.0014 prescripted in 
Ref. However, on this wav the uncertainty of the proton radius is significantly larger 
than in both Mainz papers B i (see Ref. i for more details). 

We have discussed all items of Table 1 and we present now some conclusions. The 
scattering value has an uncertainty of about 0.24 fm. We cannot present here any eventual 
figure because Wong performed his evaluation with the most important part of data but 
not with all of them. Not all corrections to the data were included. The result should be 
close to the Wong value (0.877(24) fm) but it needs more analysis. However data are quite 
old and practically it is not possible to reevaluate them. Hopefully, a new measurement 
in Mainz is going to be done. The Lamb shift examination can give a better value. We 
have to mention that the experiment claimed as the most accurate one (uncertainty is 2 
kHz) is under some criticism and the average value of all other measurements gives now 
some larger uncertainty (about 3 kHz). The uncertainty for the proton radius extracted is 
0.15 fm (the equal portions arise from the ground state Lamb shift QED theory and from 
the optical experiments) . We also would not like to give any value because it is a running 
value now: some data are still coming. The most precise value which is included into 
our analysis is due to the Garching and Paris absolute measurements. Some parts of the 



Paris data are changing now due to recalibration of the standards Izzl. The preliminary 
value from the Lamb shift is in fair agreement with the Wong value but with twice larger 
accuracy. 

The situation now, when the uncertainties from spectroscopy measurement, QED 
calculation and the scattering data are on the same level is quite challenging and we 
believe that new results are coming in all these fields. Recently a new wave of activity 
projecting of a muonic hydrogen experiment for the Lamb shift has been arisen at the 
PSI and we also hope that will give one more value for the proton radius. We also would 
like to attract attention to the Lamb shift measurement in hydrogen-like system with a 
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moderate value of the nuclear charge Z ~ 5 — 20. It seems it is necessary to encourage 
such investigations because they can lead to experim,ental estimation of the higher order 
QED terms. That is to be helpful for the hydrogen Lamb shift theory. 

This work was supported in part by the Russian state program of Fundamental 
metrology. We would hke to thank T. W. Hansch and the whole hydrogen team of the 
MPQ for the support, hospitality and stimulating discussions. 
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